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ABSTRACT 

We present a survey of 2,148 galaxy redshifts from the VLT LBG Redshift Survey 
(VLRS) , a spectroscopic survey of z ~ 3 galaxies in wide fields centred on background 
QSOs. We have used deep optical imaging to select these galaxies via the Lyman 
break technique. Spectroscopy of the Lyman Break Galaxies (LBGs) was then made 
using the VLT VIMOS instrument, giving a mean redshift of z w 2.78. To make a 
definitive LBG clustering analysis, we have combined the VLRS redshifts with the 
813 Keck LBG redshifts of Steidel et al, with the statistical power of VLRS at large 
scales complementing the accuracy of the Keck sample at small scales. From the semi- 
projected correlation function, w p (cr), for the VLRS and combined surveys, we find 
that the results are well fit with a single power law model for the real space correlation 
function, £(r), with clustering scale lengths of respectively tq = 3.32 ±0.41 and 3.75 ± 
0.24 /i _1 Mpc. We note that the corresponding combined £(r) slope is flatter than for 
local galaxies at 7 = 1.55 ± 0.09 rather than 7 = 1.8. This flat slope is confirmed by 
the z-space correlation function, £(s), and in the range 10 < s < 100 /i Mpc the 
VLRS shows an « 2.5cr excess over the KCDM linear prediction. This excess may be 
consistent with recent evidence for non-Gaussianity in clustering results at z ~ 1. We 
then analyse the LBG z-space distortions using the 2-D correlation function, £(cr, 7r), 
finding for the combined sample a large scale infall parameter of (3 = 0.32 ± 0.20 and 
a velocity dispersion of yj (tuj) = 540t2ookms~ 1 . Fixing this velocity dispersion, we fit 
£(cr, 7r) for Q m and /3 and break their degeneracy using low-z data to find fl m (z = 0) = 
0.301 

18- Finally, based on our measured /3, we are able to determine the gravitational 
growth rate, finding a value of f(z = 3) = 0.83 ± 0.46 (or fa s = 0.22 ± 0.12), which 
is the highest redshift measurement of the growth rate via galaxy clustering and is 
consistent with Einstein gravity and ACDM. 

Key words: galaxies: kinematics and dynamics - cosmology: observations - large- 
scale structure of Universe 



* Based on data obtained with the NOAO Mayall 4m Telescope at Kitt Peak National Observatory, USA (programme ID: 06A- 
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1 INTRODUCTION 

The large scale structure of matter presents a crucial guide 
in understanding the nature and evolution of the Uni- 
verse. In ACDM, structure in the Universe grows hierar- 
chica l ly through gravitational instability (e.g. iMo fc White] 
1 19961 ; Ijenkins et aflll99St ISpringel et al]|2006l ) and testing 
this model requires the measurement of the matter clus- 
tering and the growth of structure a c ross cosmic time (e.g. 
ISpringel et afll2005l : lOrsi et alj|2008l ; iKim et alj|2009h . We 
are limited however in our ability to trace the structure of 
mass given that the ACDM model suggests that « 80% of 
the mass density of the Universe is in the form of dark mat- 
ter. 

Although large photometric surveys are beginning to 
map the over all matter density distribution via its lensin g 
signature (e.g. iMassev et al. I l2007l ; iHildebrandt et al.ll2012f ) , 
at present the key tool in the statistical analysis of the dis- 
tribution of matter in the Universe remains the study of the 
clustering statistics of selected galaxy populations. A given 
galaxy population traces the peaks in the matter distribu- 
tion and hence provides a biased view of the matter density, 
which nevertheless can be used to follow the overall growth 
of structure. 

At low-redshift, magnitude limited galaxy samples have 
provided significant tools in probing the clustering prop- 
erties of the galaxy population (e.g. iNorberg et alj l2002i ; 
lHawkins et~alT 120031 ") ■ Moving to higher redshifts, photo- 
metric selections are required to isolate the required red- 
shift range, for example the Luminous Red Galaxy (LRG), 
Emission Line Galaxy (ELG), BzK, Extremeley Red Ob- 
ject (ERO), Distant Red Galaxy (DRG) selections. At 
> 2, identifying galaxy populations is primarily r eliant 



on the Lyman B r eak f e .g. ISteidel et al.lll996l . Il999l . 120031 ; 
lAdelberger et al.l 12003); iFranx et al.1 IgOOaT and the Ly-a 



emitter (LAEs; e g. ICowie fc Hulll998l ; iGawiser et al.|[200r3 . 
120071 ; lOuchi et al.ll2008l ) selections. 

Lyman break galaxies (LBGs) facilitated for the first 
time observational estimates of the st ar-formation history 
of the Universe o ut to high redshifts (|Madau et al.l 1 19961 ; 
ISteidel et alJ[l999T ) . With the possible exception of observa- 
tions made at sub-millimetre (submm) wavelengths, where 
a strong negative fc-correction results in a redshift distribu- 
tion dominated by high-z sources (e.g. Ilvison et al.ll2002l ; 
ISmail et al.1 12002] ; IChapman et al.1 120051 ). this use of broad- 
band colour criteria to preselect high-z galaxy candidates 
is the most effective tool av ailable for surveyin g the galaxy 
population at high redshift (|Steidel et al.|[2003h . 

Importantly too, w hereas submm sources are unusual 
ultra-luminous objects (|Chapman et al.l l2005h with rela- 
tively low densities, LBGs represent a much larger popula- 
tion of comparatively 'normal' galaxies. For example, their 
luminosities, star-formation rates and level of obscuration 
are found to be in line with those of l ocal star-forming galax- 
ies (e.g. ISteidel et all Il99rj. Il999l; lLowenthal et al.l Il997t 



iPettini et al.ll200ll ; iMagdis et alj|20ich . Additionally 



0133), the NOAO Blanco 4m Telescope at Cerro Tololo Inter- 
American Observatory, Chile (programme IDs: 03B-0162, 04B- 
0022) and the ESO VLT, Chile (programme IDs: 075.A-0683, 
077.A-0612, 079.A-0442). 
f E-mail: rmbielby@gmail.com (RMB) 



pared to the selection of LAEs , which offer an alternative 
selection of galaxies at z > 2 via ground-based observations, 
the LBG selection allows the identification of galaxies across 
a br oad and cont i guous range of redshifts. 

ISteidel et al.l ( 20031 ) presented a large survey of LBGs 
in the re dshift range 2.5 < z < 3. 5, identifyin g ~ 800 such 
galaxi es. Ida Angela et al.l (|2005al ) took the ISteidel et al.l 
(2003) sample and used the clustering properties of the LBG 
population to constrain the cosmological density parameter, 
f2 m , and the bulk motion properties of the large scale struc- 
ture at z ~ 3. They measured a clustering amplitude for 
the LBGs of r = 4.48tg; 1 g and based on the 2D cluster- 
ing of the galaxies place constraints on the infall param- 
eter of /3(z = 3) = 0.251o'q6 and on the mass density of 
f2m(2 = 0) = 0.551q' 16 . These results assume a velocity dis- 
persion, which is degenerate with the infall parameter, of 
400kms _1 and are based on small-scales with the transverse 
separation li mited to a < 10 fe -1 Mpc and so were some- 
what limited. iBielbv et al.l (|201ll . Paper I) added to this pre- 
senting the first data from the VLT LBG Redshift Survey 
(VLRS). Adding ~ 1,000 galaxies to the 2 « 3 sample of 
ISteidel et al.l (|2003h data, they measured the clustering and 
dynamics of the 2 « 3 LBG population. 

In this paper, we add to the previous resul ts of the 
VLT LBG Redshift S urvey (VLRS) presented in iPaper J . 
ICrighton et~al] (|201ll . Paper II) and IShanks et al.l (|201lf ). 
Here we measure the clustering of the z « 3 LBG popula- 
tion. Throughout this paper, we use a cosmology given by 
H = 70kms"\ ft m = 0.3, fi A = 0.7 and a 8 = 0.8. In addi- 
tion distances are quoted in comoving coordinates in units 
of /i -1 Mpc unless otherwise stated. 



2 OBSERVATIONS 
2.1 Survey overview 

The VLRS presently encompasses nine separate fields cen- 
tred on bright z ~ 3 QSOs. The imagi ng and s pectroscopy 
for the first five fields were presented in lPaper il . The obser- 
vations described herein form the second part of the VLT 
LBG Survey. In order to facilitate an investigation of how 
z m 3 galaxies interact with gas in the intergalactic medium 
(IGM), the survey comprises observations of several target 
fields centred on bright z > 3 quasars, since features in the 
QS O spectra can provide information on the local IGM. 

IPaper J presented the first 5 fields of the survey, cen- 
tred on the following quasars: Q0042-2627 (z = 3.29), 
J0124+0044 (z = 3.84), HE0940-1050 (z = 3.05), 
J1201+0116 (z = 3.23) and PKS2126-158 (z = 3.28), here- 
after referred to by only the right-ascension component of 
these names. A spectroscopic survey of each of these quasar 
fields was carried out with the Visible Multi-Object Spec- 
trograph (VIMOS) on the European Southern Observatory's 
Very Large Telescope (VLT) in Chile (during the ESO peri- 
ods 75-79). Each field consisted of four sub-fields (individual 
pointings with the VLT spectrograph), except for HE0940 
where only three sub- fields were available at the time of their 
publication. A VIMOS pointing has a field of view of 16' x 18' 
(see H2.4.ip . therefore each quasar field covered « 32' x 36', 
or « 0.32 deg 2 , except for HE0940 which with 3 sub-fields 
covered « 0.24 deg 2 . 
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Table 1. A summary of the fields making up our z rs 3 LBG survey. The table gives the name, coordinates and redshift of the QSO on 
which the fields are roughly centred, as well as the number of subfields (individual VLT VIMOS pointings) with spectroscopic data. The 
first block of fields were presented by I Paper j , the second block are presented in this paper. 



Field 


RA a 


Dec" 


z b 


Subfields 


Reference 


Q0042-2627 


00:44:33.9 


-26:11:21 


3.29 


4 


PaDer I 


J0124+0044 


01:24:03.8 


+00:44:33 


3.84 


4 


Paper 1 


HE0940-1050 


09:42:53.4 


-11:04:25 


3.05 


3 


Paper I 


J1201+0116 


12:01:44.4 


+01:16:12 


3.23 


•1 


Paper I 


PKS2126-158 


21:29:12.2 


-15:38:41 


3.28 


1 


Paper I 


19 


Q2359+0653 


00:01:40.6 


+07:09:54 


3.23 


4 


This work 


Q0302-0035 


03:03:41.0 


-00:23:22 


3.23 


4 


This work 


Q2231+0015 


22:34:09.0 


+00:00:02 


3.02 


3 


This work 


HE0940-1050 


09:42:53.4 


-11:04:25 


3.05 


6 


This work 


Q2348-011 


23:50:57.9 


-00:52:10 


3.02 


9 


This work 


26 



a J2000 coordinates of QSO; not necessarily the exact centre of the observed field. 
b redshift of the central quasar 



Building on this initial dataset, we present the continu- 
ation of these observations since incorporating ESO periods 
81 and 82. We have added a further 6 sub-fields to HE0940, 
tripling its previous area, as well as observations of 4 new 
fields, around the quasars Q2359+0653 (z = 3.23), Q0302- 
0035 = 3.23), Q2231+0015 (z = 3.02) and Q2348-011 
(z = 3.02), with 4, 4, 3 and 9 sub-fields respectively. Ta- 
ble [T] summarises a ll the fie lds of the survey. This includes 
those presented by |Paper"H . covering 1.52 deg 2 , and those 
presented here, which take the total observed area to 3.6 
deg 2 , more than doubling the previous size. 

2.2 Imaging 

2.2.1 Observations and data reduction 

The selection of z ~ 3 LBG candidates was performed us- 
ing photometry from optical broadband imaging. The imag- 
ing data for Q2359 and Q0302 were acquired with the Mo- 
saic wide-field imager on the 4m Mayall telescope at Kitt 
Peak National Observatory (KPNO) in September 2005. 
The Q2231 data are from the Wide Field Camera on the 
2.5m Isaac Newton Telescope (INT) on La Palma, and were 
observed in August 2005. All of these observations were car- 
ried out in the U, B and R bands. 

The HE0940 and Q2348 data were acquired with the 
MegaCam imager on the 3.6m Canada-France-Hawai'i Tele- 
scope (CFHT). HE0940 was observed in the u, g, r, i and z 
bands in April 2004 as part of the observing run 2004AF02 
(PI: P. Petitjean), whilst Q2348 was observed in the u, g, r 
and z bands over the period August-December 2004 as part 
of the observing run 2004BF03 (PI: P. Petitjean). Table [2] 
gives full details of all the imaging data. The instruments, 
filters and data reduction processes are described below. 

The Mosaic imager at KPNO consists of 8 2kx4k 
CCDs arranged into an 8kx8k square. With a plate scale of 
0.26"/pixel, this gives a field of view of 36'x36'. There are 
0.5-0. 7mm gaps between the chips, corresponding to gaps 



of 9-13" on-sky, so a dithering pattern was used during the 
observations to provide complete field coverage. U, Harris 
B and Harris R filters were used, centred at 3552A, 4298A 
and 6380A respectively. 

The Mosaic image data were reduced using the mscred 
package in iraf. The reduction process is described by 
IPaper J . however we briefly outline the procedure here. Ini- 
tially a master bias frame is produced for each night's ob- 
serving. The dome flats and sky flats were then processed 
using the cedproc and mcspupil routines, subtracting the 
bias and eliminating the faint 2600-pixel pupil image arte- 
fact. The object frames were processed similarly, subtracting 
the bias and pupil image, and then were flat-fielded using the 
dome and sky flats. Bad pixels and cosmic rays were masked 
out of the science frames using the crreject, crplusbpmask 
and fixpix procedures . Finally, the SWARP software pack- 
age (|Bertin et all 12002) was used to resample and co-add 
the frames, producing a final science image. 

The Wide Field Camera (WFCam) on the INT com- 
prises 4 2kx4k CCDs. These are arranged into a 6kx6k 
block with a 2k x 2k square missing. With « 1' gaps between 
chips and a pixel scale of 0.33"/pixel, WFCam has a total 
FoV of « 34'x34' (0.32 deg 2 ); however, accounting for the 
incomplete coverage of the field, the total observing area is 
reduced to 0.28 deg 2 . 

The WFCam observations of Q2231 were made using 
the RGO U, Harris B and Harris R filters. The RGO U filter 
has a central wavelength of 358lA and a FWHM of 638A, 
making it very similar to the U band filter used at KPNO 
(centre 3552A, FWHM 631 A). The B and R band filters 
were the same as at KPNO. Therefore, given that the filters 
are so similar, we will use the same UBR selection criteria 
when identifying LBG candidates in either the Mosaic or 
WFCam datasets. 

Initial data reduction, including bias removal, flat- 
fielding and photometric calibration, was performed by the 
Cambridge Astronomical Survey Unit (CASU). Astrometry 
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Table 2. Details of imaging observations for the LBG target fields presented in this paper. 



Field 


RA 


Dec 


Instrument 


Band 


Exposure 


Seeing 


Completeness 


Dates 




(J2000) 






(ks) 




(50% Ext/PS) 




Q2359 


00:01:44.85 


+07:11:56.0 


Mosaic (KPNO) 


U 


19.2 


1.46" 


24.76/25.18 


29-30 Sep 2005 










B 


7.2 


1.45" 


25.28/25.73 












R 


6.0 


1.15" 


24.74/25.20 




Q0301 


03:03:45.27 


-00:21:34.2 


Mosaic (KPNO) 


U 


19.2 


1.34" 


24.93/25.34 


29-30 Sep 2005 










B 


6.4 


1.28" 


25.51/26.04 












R 


4.8 


i i nil 
1.19 


24.59/25.17 




Q2231 


22:34:28.00 


+00:00:02.0 


WFCam (INT) 


U 


54.0 


1.23" 


25.08/25.52 


30 Aug 2005 










13 


13.2 


1.01" 


25.88/26.12 












R 


19.2 


1.01" 


24.75/25.24 




HE0940 


09:42:53.06 


-11:02:56.9 


MegaCam (CFHT) 


u 


6.8 


0.99" 


25.39/25.93 


14, 21-27 Apr 2004 










9 


3.1 


0.86" 


25.54/26.05 












r 


3.7 


0.85" 


25.08/25.65 




Q2348 


23:50:57.90 


-00:52:09.9 


MegaCam (CFHT) 


u 


9.9 


0.78" 


25.97/26.62 


19-20 Aug, 










9 


5.5 


0.79" 


25.71/26.29 


7-10 Nov, 










r 


1.1 


0.75" 


25.22/25.80 


15 Dec 2004 



calibration and exposure stacking was performed using the 
SCAMP and SWARP packages. 

MegaCam is a considerably larger instrument than Mo- 
saic or WFCam, with 36 2k x 4.6k CCDs arranged into a 
18kxl8k grid. With a high resolution of 0.187" /pixel this 
gives a 1° x 1° FoV. The observations of HE0940 and Q2348 
were carried out using the u*, g' and r' filters (referred to 
in this paper as u, g and r for convenience), with respec- 
tive central wavelengths of 3743A, 4872A and 6282A. These 
filters are significantly different from the filters used in our 
other fields and alternative LBG selection criteria are re- 
quired for the MegaCam datasets. 

2.2.2 Photometry 

Photometric zeropoints for the imaging fields were deter- 
mined from standard star observations carried out as part 
of each of the imaging runs. The standard star fields were 
reduced in the same way as the science frames to ensure 
consistency. Source detection in the science images was per- 
formed with SExtractor (Bertin & Arnouts, 1996), using a 
1.5a detection threshold and a 5-pixel minimum size. 

The Uvega, Bvega and Rvega band galaxy number 
counts in the Q0301 (diamonds), Q2231 (triangles) and 
Q2359 (squares) LBG fields are shown in the top panels of 
Fig. [1] Stars were removed from these counts at magnitudes 
brighter than « 22 using a limit on the measured half-light 
radius of the sources. At fainter magnitudes, no attempt to 
remove stars from the counts was made, as the smallest ex- 
tended sources become unresolved at the PSF of our fields 
at such magnitudes. We also show completeness estimates 
for each image in each field. These are estimated by plac- 
ing simulated sources at random positions in a given image 
and measuring the fraction that are successfully extracted 
with SExtractor (using the same extraction parameters as 
used to create the full catalogues). In each case we estimate 
the completeness using both simulated point-sources and ex- 
tended sources, where the extended sources are modelled 



by a de Vaucouleurs r 1//4 profile with a half-light radius of 
r i/2 = 0.3". In both cases, the simulated source is convolved 
with the image PSF before being added to the observation. 

The results of the completeness estimates for the 
U BRvega filter fields are shown in the lower panels of Fig. [T] 
The same symbols as the top panels are used for the differ- 
ent fields, whilst the dashed curves show the completeness 
estimates based on the extended sources and the solid curves 
show the completeness for the simulated point-sources. The 
50% limits completion estimates are given in Table [21 Com- 
paring the completeness measurements across the fields, the 
measurements are relatively consistent with the imaging in 
each field reaching comparable depths. 

We show the galaxy number counts (top panels) and 
completeness estimates (lower panels) for the MegaCAM 
fields in Fig. [2] Again the symbols are consistent between 
top and lower panels with the diamonds showing the results 
for the HE0940 field and the triangles showing the Q2348 
field. As before, the solid lines in the lower panels show the 
completeness estimates for the point-like sources and the 
dashed lines show the same for the extended sources (which 
use the same de Vaucouleurs profile as used for the U BRvega 
fields). Comparing the two fields to each other, the depths 
reached are comparable in each band, although the HE0940 
is marginally less deep in the u band by « 0.5 mag. 

2.3 Candidate selection 

2.3.1 UBR selection 

In the Q2359, Q0302 and Q2231 fields, we selected LBG 
candidates based on their U, B and R p hotometr y. The cri- 
teria used were th e same as thos e used by |Paper"H . which are 
based on those of ISteidel etall (|2003l L There are 4 groups 
to the selection, designated lbg_pri1, lbg_pri2, lbg_pri3 
and lbg_drop and defined as follows: 

LBG_PRIl 
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Figure 1. Upper panel: Galaxy number counts in the U (left), B (middle) and R (right) band imaging for the fields Q0301 (diamonds), 
Q2231 (triangles) and Q2359 (squares). Lower panel: Estimated completeness for each of the above bands based on simulated point 
sources (solid lines) and Vaucouleurs profile sources (dashed lines). 



• 23 < R < 25.5 

• 0.5 <(U-B) <4.0 

• (B — R) < 0.8(1/ - B) + 0.6 

• (B - R) < 2.2 

LBGJPRI2 

• 23 < R < 25.5 

• {U-B)> 0.0 

• (B — R) < 0.8(U - B) + 0.8 

• -1< (B-R) < 2.7 

• ^ LBG_PRIl 
LBG_PRI3 

• 23 < R < 25.5 

• -0.5 <{U-B) <0.0 

• -1.0 < (B — R) < 0.8(U — B) + 0.6 

• ^ {LBG_PRI1,LBG_PRI2} 
LBG_DROP 

• 23 < R < 25.5 

• 0.5 <(B-R) < 2.2 

• No detection in U 

The first 3 groups represent an order of priority - 
that is, LBG-PRll candidates are considered more likely to 
be z « 3 LBGs than e.g. LBG_PRI3 candidates. This is be- 
cause whereas LBG_PRll tends to select outliers in the UBR 
colour-colour plot, the lower priority groups select objects 
increasingly close to the colour region populated by stars and 



lower redshift galaxies, and therefore suffer from increased 
contamination from lower-redshift interlopers. 

The fourth group is somewhat separate, being for galax- 
ies which are not detected in the U band. Such sources may 
be excellent LBG candidates, since it may be that the pres- 
ence of the Lyman limit in the U band has made the galaxy 
extremely faint in this band, such that it 'drops out' below 
the magnitude limit. However, the lbg_drop population is 
also likely to suffer from contamination, in this case because 
objects with no counterpart in 1 of the 3 bands have a higher 
chance of being spurious sources. 

Figs. O 2] and show UBR colour-colour plots for 
Q2359, Q0302 and Q2231, respectively. In each plot, the 
LBG_PRll, LBG_PRI2, LBG_PRI3 and lbg_drop candidates are 
indicated. A model colour-redshift track is also plotted, 
showing the expected evolutionary path of a star-forming 
galaxy from z = to z = 4. This was produced using the 
Bruzual & Chariot (2003) model, assuming a Chabrier IMF 
and an exponential SFR with e-folding time r = 9 Gyr. The 
model indicates that our selection criteria (across all the 
priority groups) is predicted to isolate galaxies in the range 
fn 2.5 < z < 3.8. It also suggests that, of the sources that 
are confirmed as high-redshift LBGs, the LBG_PRl3s should 
typically be at a lower redshift than the LBG_PRl2s, which in 
turn should be at lower redshift than the LBG-PRlls. |Paper"ll 
noted that this trend was detected in their LBG sample. 
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Figure 2. Number counts (top) and completeness estimates (bottom) in the u (left), g (middle) and r (right) bands from the McgaCAM 
imaging on the HE0940-1050 and Q2348-011 fields. 




Figure 3. UBR colour-colour plot showing candidate selection Figure 4. As for Fig. [3] but for Q0302. 

in Q2359. Objects selected as LBG_PRll, LBG_PRl2, LBG_PRl3 and 
LBGJDROP candidates are shown in different colours as indicated 
by the legend. The LBG_DROP candidates have been placed at 
U — B = 4.5. The contours show the colour distribution of the 

rest of the objects in the field. 2.3.2 ugr selection 

In HE0940 and Q2348, LBG candidates were selected based 
on ugr photometry. We have therefore adapted the cri- 
teria outlined above to account for the different colour 
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LBGPRH 
LBG_PRI2 
X LBGPRI3 
LBG DROP 





LBGPRH X 
LBG_PRI2 
LBG PRI3 X 
LBG DROP 



7 



-2-10123 
g-r 



Figure 5. As for Fig. El but for Q2231. 



Figure 7. As for Fig. El but for Q2348. 



LBGPRH 
LBG_PRI2 
LBG PRI3 X 
LBG DROP 




Figure 6. A ugr colour-colour plot showing the selection of can- 
didates in HE0940. Symbols are as in Fig. El 



bands. Again candidates were selected as either lbg_pri1, 
lbg_pri2, LBG_PRI3 or lbg_drop, defined as follows: 

LBG.PRIl 

• 23 < r < 25 

• 1.4 <(u-g)< 4.0 

• -0.36 < (g-r) < 0.96 
. (g-r) < (u- g) -1.88 

LBG-PRI2 



• (g - r) < (u - g) - 1.44 

• ^ LBG_PRIl 

LBG-PRI3 

• 23 < r < 25 

• (g-r) < (u-g) -0.7 

• -0.45 < (g - r) < 0.2(u - g) + 0.1 

• ^ {LBG_PRI1,LBG_PRI2} 
LBG_DROP 

• 23 < r < 25 

• -0.36 < (g - r) < 0.96 

• No detection in u 



• 23 < r < 25 
. (u-g) > 1.0 

• -0.36 < (g-r) < 0.96 



Figs. [B] and [7] show the resulting candidates on ugr 
colour-colour plots, and Table [3] gives the numbers and sky 
densities of candidates in all 5 LBG fields. 



2.4 Spectroscopy 

2-4-1 Observations 

The LBG candidates were targeted in spectroscopic follow- 
up observations with the VLT VIMOS spectrograph be- 
tween September 2008 and December 2009, with programme 
IDs081.A-0418(B) (Q2231), 081.A-0418(D) (Q2359), 081.A- 
0418(F) (Q0302), 082.A-0494(B) (HE0940) and 082. A- 
0494(D) (Q2348). The observations were done during dark 
time in generally good conditions with a typical seeing of 
w 1" and an airmass of 1.0-1.3. Table [4] gives details of all 
the fields observed. 

The VIMOS instrument |Le Fevre et all 120031 ') com- 
prises four separate CCDs, each with a field of view of 7' x 8'. 
These four arms are arranged in a 2 x 2 grid with a ~ 2" gap 
between each CCD, giving a total 16' x 18' FoV as quoted 
previously. Of this 288 arcmin 2 field, 224 arcmin 2 is covered 
by the detector. 

Our observations utilised the low resolution blue 
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Table 3. Numbers and sky densities of z 3 LBG candidates in each of the fields. For each priority class, the first column shows the 
total number of candidates selected and the second, italicised column shows the density in arcmin™ 2 . The last two columns show the 
figures for all z « 3 LBG candidates. 



Field LBG_PRll LBG_PRl2 LBG_PRl3 LBGJDROP All 



Q2359 


795 


0.61 


1,130 


0.87 


549 


0.42 


2,235 


1.72 


4,709 


3.63 


Q0302 


891 


0.69 


1,227 


0.95 


433 


0.33 


1,014 


0.78 


3,565 


2.75 


Q2231 


748 


0.65 


948 


0.82 


424 


0.37 


514 


0.45 


2,634 


2.29 


HE0940 


2,808 


0.78 


1,494 


0.42 


4,347 


1.21 


6,146 


1.71 


14,795 


4.11 


Q2348 


1,843 


0.51 


1,624 


0.45 


4,808 


1.34 


1,850 


0.51 


14,795 


2.81 



(LR_Blue) grism and the order sorting blue (OS _Blue) filter, 
resulting in a wavelength range of 3700A-6700A, blazing at 
~4000A. This wavelength range is ideal for our survey, de- 
tecting the Lya line at 2.0 < 2 < 4.5. The resolving power 
of the spectrograph in this configuration is R = 180 assum- 
ing a 1" slit (as used in these observations), which gives a 
resolution element of AA « 22A at the blaze wavelength. 
The spectral dispersion is 5.3A/pixel. 

The slit masks were designed using the vmmps software 
which is standard for VIMOS observations. The aims for 
mask design are (a) to maximise the number of observed tar- 
gets, (b) to favour higher-priority targets and (c) to ensure 
slits are of sufficient size to allow a robust sky subtraction. 
Since these aims are frequently in conflict with one another, 
the mask design process is one of attempting to optimise the 
observations to satisfy all three as well as possible. Point (c) 
is addressed by setting a minimum slit length of 8" (40 pix- 
els given the pixel scale of 0.205" /pixel)- Slit length was 
increased as much as possible where such an increase would 
not prevent the observation of an additional target — that 
is, where it did not conflict with point (a). 

With the LR_Blue grism, each spectrum spans 640 pix- 
els along the dispersion axis. Assuming a 40 pixel slit width 
as specified above, this would allow for a possible total of 
over 300 slits on the full 4kx2k detector. This is however not 
practically achievable given the density of LBG candidates, 
and is hampered further by the need to select high-priority 
candidates (point b), which have an even lower sky density. 
Our final slit masks therefore typically contain some 50-70 
slits per quadrant. 



2.4.2 Data reduction 

The spectroscopic data have all been reduced using the VI- 
MOS ESOREX reduction pipeline. Using bias frames, flat 
fields and arc lamp exposures taken for each mask during 
each observing run, the pipeline generates bias-subtracted, 
flat-fielded, wavelength-calibrated science frames consist- 
ing of a series of 2D spectra. Following iPaper j we use 
the imcombine procedure in iraf to combine the reduced 
frames from each observing block, generating a master sci- 
ence frame for each quadrant of each field. When combin- 
ing the frames we use the crreject mode, designed to re- 
move cosmic rays by rejecting pixels with significant posi- 
tive spikes. We have also used the avsigclip rejection mode 
with a rejection threshold of 3a, and find that our results 
are not significantly affected, suggesting that our results do 



not depend strongly on the parameters used to combine the 
science frames at this stage. 

We extract ID spectra from the reduced, combined 2D 
spectra using the idl routine specplot. One-dimensional 
object and sky spectra are found by averaging across the re- 
spective apertures, and the sky spectrum is then subtracted 
from the object spectrum to give a final spectrum for the 
object. 

In some cases there remain significant sources of con- 
tamination in the final object spectrum. These can arise 
from bad pixels, either in the object or sky aperture, from 
contamination from the zeroth order from other slits, or 
more frequently from the bright sky emission lines [0 1] 
5577A, [Nal] 5990A and [Ol] 6300A; in either case, the re- 
sulting contamination may manifest itself as either a posi- 
tive or a negative spike in the spectrum. Such artefacts are, 
however, easily spotted during a routine inspection of the 
2D spectrum. 

2.5 Identification of targets 

Every source targeted for spectroscopic observation is in- 
spected visually, in both the 2D and ID spectra, to deter- 
mine where possible a redshift and classification. Sources 
are classified as either z m 3 Lyman-break galaxies, low- 
redshift galaxies, QSOs or Galactic stars. The LBGs are 
divided into those showing Lya emission (designated LBe) 
and those showing Lya absorption (LBa). QSOs are deter- 
mined by the presence of broad emission features, in par- 
ticular Lya and CIV. Stars are classified by comparison to 
template spectra: in particular we check for A, F, G, K and 
M stars. 

In determining the redshift and classification the spec- 
tral feature primarily used in the case of LBGs is the Lya 
emission/absorption line at 1216A; for lower redshift galax- 
ies it is the [0 11] emission line at 3727A. In addition to these, 
some of the following features are used: 
For z « 3 LBGs: 

• Lyman limit, 912 A; 

• Ly/3 emission/absorption, 1026A 

• O VI 1032A, 1038A; 

• Lya forest, <1215.67A; 

• Lya emission/absorption, 1215. 67A; 

• Inter-stellar medium (ISM) absorption lines: 

- Sill 1260.4A; 

- Ol+Sill 1303A; 

- CII 1334A; 
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Table 4. Details of spectroscopic observations for the LBG target fields presented in this paper. 



Field 


Subficld 


RA a 


Dec a 


Exposure 6 


Airmass 


Seeing 


Dates 


Q2359 


fl 


00:01:09.94 


+07:03:26.8 


10 


1.1 - 


1.2 


0.6 - 


1.7" 


23-25 Sep 2008 


Q2359 


f'2 


00:01:12.92 


+07:16:39.2 


10 


1.1 - 


1.3 


0.6 - 


1.3" 


3, 20-21 Oct 2008 


Q2359 


f3 


00:02:11.50 


+07:15:33.9 


10 


1.2 - 


1.4 


0.6 - 


1.3" 


3, 21-25 Nov 2008 


Q2359 


fl 


00:02:12.89 


+07:02:19.1 


10 


1.2 - 


1.3 


0.5 - 


1.1" 


26-30 Nov 2008 


Q2231 


fl 


22:34:28.19 


-00:06:03.3 


10 


1.1 - 


1.2 


0.5 - 


1.0" 


23, 28 Oct 2008 


Q2231 


f2 


22:34:28.55 


+00:06:13.2 


10 


1.1 - 


1.2 


0.4- 


0.9" 


21-22 Oct 2008 


Q2231 


13 


22:33:39.51 


-00:06:10.8 


10 


1.1 - 


f .2 


0.5 - 


1.0" 


3 Aug; 27 Jul 2008 


Q0302 


fl 


03:04:20.12 


-00:14:28.8 


12 


1.1 — 


f .3 


0.7 — 


1.2" 


23, 31 Oct 2008 


Q0302 


f2 


03:03:10.27 


-00:16:18.7 


12 


1.1 — 


f .2 


0.7 — 


1.5" 


21-23 Nov 2008 


Q0302 


13 


03:03:15.41 


-00:30:40.0 


12 


1.1 — 


f .2 


0.7 — 


1.5" 


25-26 Nov 2008 


Q0302 


fl 


03:04:15.56 


-00:28:59.1 


12 


1.1 — 


1,2 


0.7 — 


1.2" 


24 Sep; 1, 7 Oct 2008 


HE0940 


f4 


09:42:10.00 


-10:54:30.3 


11.2 


1.0 - 


1.3 


0.8 - 


1.5" 


1 Feb 2009 


HE0940 


f5 


09:43:07.47 


-11:24:50.3 


11.2 


1.0 - 


1.3 


0.7- 


1.4" 


3 Feb 2009 


HE0940 


ffi 


09:41:59.99 


-11:24:50.4 


11.2 


1.0 — 


f .2 


0.5 — 


1.2 


20-21 Feb 2009 


HE0940 


17 


09:44:14.99 


-11:24:49.9 


11.2 


1.0 - 


1.2 


0.5 - 


1.3" 


22, 24 Feb 2009 


HE0940 


i'S 


09:43:21.49 


-10:41:00.5 


11.2 


1.0 - 


1.2 


0.5 - 


1.0" 


26-27 Feb 2009 


HE0940 


ffl 


09:42:09.99 


-10:40:59.8 


11.2 


1.0 - 


1.3 


0.5 - 


1.2" 


2 Feb 2009 


Q2348 


fl 


23:51:50.08 


-00:54:21.9 


11.5 


1.1 - 


1.2 


0.5 - 


1.7" 


23-25 Jul 2009 


Q2348 


12 


23:50:45.09 


-00:54:22.2 


11.5 


1.0 - 


l.f 


0.5 - 


1.0" 


19-20 Jul 2009 


Q2348 


13 


23:49:40.07 


-00:54:22.6 


11.5 


1.0 - 


1.2 


0.4- 


0.8" 


27 Jul 2009 


Q2348 


ff 


23:51:50.12 


-00:37:31.6 


11.5 


1.1 - 


1.2 


0.5 - 


1.5" 


20-21 Aug 2009 


Q2348 


f'5 


23:50:45.05 


-00:37:31.5 


11.5 


1.0 - 


1.2 


0.7- 


1.4" 


16-20 Sep 2009 


Q2348 


ffi 


23:49:40.00 


-00:37:32.0 


11.5 


1.1 - 


1.2 


0.8 - 


1.3" 


24-25 Sep 2009 


Q2348 


f7 


23:51:50.12 


-01:07:31.4 


11.5 


1.1 - 


1.2 


0.7- 


1.0" 


12, 20 Oct 2009 


Q2348 


1'8 


23:50:45.00 


-01:07:32.0 


11.5 


1.1 - 


1.3 


0.7- 


1.3" 


22 Nov, 10 Dec 2009 


Q2348 


f9 


23:49:40.00 


-01:07:32.0 


11.5 


1.1 - 


1.3 


0.8 - 


1.5" 


15-22 Nov 2009 



a J2000 coordinates of subficld centre 
b in ks 



the quality parameters indicate that the redshift is based on 
the following features: 

• Q = 0.5 — a spectral break with some weak Lya emis- 
sion/absorption and low-SNR ISM absorption features, or 
strong Lya emission but with no detected continuum 

• Q — 0.6 — a spectral break with high-SNR Lya emis- 
sion/absorption plus low-SNR ISM absorption features 

• Q — 0.7 — a spectral break with high-SNR Lya emis- 
sion/absorption plus unambiguous, high-SNR ISM absorp- 
tion features 

• Q ^ 0.8 — a spectral break with high-SNR Lya 
emission/absorption plus high-SNR absorption and emission 
lines 



2.6 LBG sample 

2.6.1 Sky densities and redshift distributions 

In total, the VLRS now consists of 2,148 LBGs in 10,080 
arcmin 2 (a density of 0.21 arc min -2 ). We show the sky dis- 
tribution of LBGs from both IPaper J . (open grey circles) 
and this paper (filled black circles) for all nine VLRS fields 
in Fig. [8] Known QSOs in the fields are also plotted (cyan 
stars). The total numbers of Q ^ 0.5 sources identified in 
each of the 5 fields presented here are given in Table [lOl We 
present examples of the first six LBGs in each of the fields in 



- Si IV doublet 1393A & 1403A; 

- Sill 1527A; 

- Fell 1608A; 

- A1II 1670A; 

• C IV doublet absorption, 1548-1550A. 
For low-z galaxies: 

• CN absorption 3833A; 

• K-band absorption 3934A; 

• HK break 4000 A; 

• emission 4102A; 

• H/3 emission/absorption 486lA; 

• O ill emission 4959A; 

• O III emission 5007A; 

The presence of the HK break causes these interlop- 
ers to appear fairly frequently in our spectroscopic samples, 
since these features mimic the Lyman break on which our 
selection is based. The ISM absorption features listed above 
for LBGs are therefore of considerable importance in iden- 
tifying genuine z ~ 3 galaxies. For every target which is 
identified, we assign a quality parameter to the redshift de- 
termined, in the range ^ Q ^ 1. A quality of Q ^ 0.4 in- 
dicates that a possible redshift has been determined, but is 
not considered a robust measurement. Above this, for LBGs, 
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Table 5. Example LBG identifications in the Q2359+0653. 



ID 


R.A. 


Dec. 


U -B 


B — R 


R 


^Lya. 


ZlSM 


QlD 




(J2000) 














VLRS J000139.85+070221.66 


0.4160563 


7.0393505 


0.63 


0.76 


23.5500 


2.4762 


2.4682 


0.5 


VLRS J000133.54+070127.57 


0.3897395 


7.0243263 


0.13 


0.15 


25.3600 


2.5707 


2.5603 


0.5 


VLRS J000118.84+070106.55 


0.3285175 


7.0184855 


1.29 


1.21 


23.8900 


3.0374 


3.0294 


0.5 


VLRS J000131. 05+070106.56 


0.3793770 


7.0184898 


0.58 


0.65 


25.3900 


2.7910 


2.7967 


0.5 


VLRS J000141. 27+070106.35 


0.4219468 


7.0184293 


1.51 


0.12 


24.3800 


2.6508 


2.6428 


1.0 


VLRS J000140.69+070044.11 


0.4195270 


7.0122533 


0.47 


0.65 


23.7800 


2.8162 


2.8082 


0.5 



Table 6. Example LBG identifications in the Q0301-0035. 



ID 


R.A. 


Dec. 


U - B 


B - R. 


R 


^Lyct 


ZISM 


QlD 




(J2000) 














VLRS J030434.85-001549.27 


46.1452103 


-0.2636854 


0.59 


0.61 


24.4700 


2.6132 


2.6041 


0.7 


VLRS J030435.40-001607.15 


46.1474953 


-0.2686527 


1.56 


1.45 


23.4100 


2.5969 


2.6157 


0.7 


VLRS J030439.49-001619.35 


46.1645317 


-0.2720422 


1.08 


0.62 


24.5100 


2.9570 


2.9490 


0.5 


VLRS J030438.22-001647.63 


46.1592560 


-0.2798966 


-0.22 


0.37 


23.8600 


2.7098 


2.7292 


0.6 


VLRS J030435.86-001654.14 


46.1494179 


-0.2817046 


0.80 


0.92 


25.0700 


2.8887 


2.8807 


1.0 


VLRS J030426. 38-001701. 38 


46.1099281 


-0.2837157 


0.21 


0.56 


24.4900 


2.4651 


2.4571 


0.6 



Table 10. Numbers of objects in each target field spectroscopi- 
cally identified (at Q ^ 0.5) as either high-z LBGs, z < 2 galaxies, 
QSOs or stars. 



Field 


2^3 LBGs 


2 < 2 galaxies 


QSOs 


Stars 


Q2359 


180 (0.20 arcrnin" 2 ) 


127 


3 


8 


Q0302 


204 (0.23 arcmin" 2 ) 


131 


10 


13 


Q2231 


114 (0.17 arcmin- 2 ) 


128 


6 


18 


HE0940 


373 (0.28 arcmin" 2 ) 


260 


3 


18 


Q2348 


335 (0.17 arcmin" 2 ) 


183 


8 


48 


Total 


1,206 (0.21 arcmin" 2 ) 


829 


30 


135 



Table 11. Rcdshift distribution statistics for spectroscopically 
confirmed, Q > 0.5 LBGs in the 5 observed fields. In each case 
the mean redshift z (with standard error), median redshift z and 
standard deviation a of the distribution is given. 



Field 


z 


z 


(7 


Q2359 


2.77 + 0.02 


2.72 


0.33 


Q0302 


2.64 + 0.02 


2.61 


0.28 


Q2231 


2.65 + 0.03 


2.62 


0.29 


HE0940 


2.81 + 0.02 


2.77 


0.33 


Q2348 


2.87 + 0.02 


2.87 


0.36 


Total 


2.78 + 0.01 


2.74 


0.34 



Tables [5] to [9] The full tables will be made available online 
at Ihttp: / /star- www.dur.ac.uk/~bielby /vlrs / 1 

For the z > 2 sample, the VLRS now consists of 1,016 
Q = 0.5, 545 Q = 0.6, 334 Q = 0.7, 153 Q = 0.8 and 100 
Q = 0.9 galaxies. 

Fig. [9] shows the n(z) distributions of all sources with 
measured redshifts in each of the 5 LBG fields. The figure 
shows that LBGs in HE0940 and Q2348, where LBGs were 
selected in ugr, have higher average redshifts than in the 
U B R-selected fields, suggesting that the ugr criteria bias 
the selection toward higher z. It is also notable from Table 
[10] that the ugr selection appears to include more Galactic 
stars. Future ugrr-selected LBG surveys may wish to alter 
our colour criteria to better avoid stellar interlopers. 

Fig. [9] also shows n(z) for the subsets of sources with 
Q ^ 0.6 and Q 0.7. We note that in any given field, the 
distributions of sources at Q ^ 0.5, Q ^ 0.6 or Q ^ 0.7 
are approximately the same — the LBGs with higher ID 
qualities are not skewed to lower or higher redshift, for ex- 



ample — suggesting that the redshift distributions shown 
are fairly robust. The average redshifts and standard devi- 
ations are given in Tab. 1111 The redshift distribution of the 
full LBG sample has a mean redshift of 2.78 ± 0.01 and a 
standard deviation of 0.34, and is shown in Fig. 1101 

Fig. [TTJ shows, as anticipated in i )2.3l that candi- 
dates selected as LBG_PRll lie at higher redshift than the 
LBG_PRI2 candidates, which are in turn at higher redshift 
than LBG_PRl3s. Quantitatively, we find that the LBG_PRlls 
have a mean redshift of z = 2.83 ± 0.02, the lbg_pri2s have 
z = 2.71 ± 0.02 and the lbg_pri3s have z = 2.61 ± 0.02. 
The lbg_drop candidates are shown in a separate panel for 
clarity, and have the highest mean redshift of all the groups, 
with z = 2.93 ± 0.02. 
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Table 7. Example LBG identifications in the HE0940-1050. 



ID 


R.A. 


Dec. 


u-g 


g-r 


r 


^Lya. 


Z ISM 


QlD 




(J2000) 














VLRS J094225.83-105744.50 


145.6076355 


-10.9623623 


2.57 


0.35 


23.6400 


2.8804 


2.8810 


0.5 


VLRS J094240.69-105753.44 


145.6695251 


-10.9648447 


1.88 


0.35 


24.1100 


3.1456 


3.1376 


0.5 


VLRS J094220.01-105900.05 


145.5833588 


-10.9833469 




-0.14 


24.3800 


2.2010 


2.1930 


0.5 


VLRS J094217.39-105923.95 


145.5724640 


-10.9899855 




0.16 


23.8500 


2.5153 


2.5073 


0.5 


VLRS J094217.51-105935.92 


145.5729675 


-10.9933100 




0.79 


24.2600 


2.8139 


2.8119 


0.6 


VLRS J094242. 29-110121. 16 


145.6762085 


-11.0225439 


0.74 


-0.03 


23.9900 


2.4613 


2.4595 


0.5 



Table 8. Example LBG identifications in the Q2231+0015. 



ID 


R.A. 


Dec. 


U - B 


B- R 


R 




ZISM 


QlD 




(J2000) 














VLRS J223439.00+000341.29 


338.6625061 


0.0614693 




0.72 


24.1500 


2.8428 


2.8291 


0.6 


VLRS J223459.87+000308.07 


338.7494507 


0.0522424 


1.01 


0.78 


23.7400 


2.4879 


2.4789 


0.5 


VLRS J223450.20+000232.38 


338.7091675 


0.0423284 




1.92 


23.7800 


2.8934 


2.8927 


0.7 


VLRS J223459.03+000051.79 


338.7459717 


0.0143855 




0.54 


24.3200 


2.8037 


2.7957 


0.9 


VLRS J223442. 76-000028.59 


338.6781616 


-0.0079414 


0.32 


1.01 


23.6900 


2.1897 


2.1817 


0.5 


VLRS J223447.81-000041.63 


338.6991882 


-0.0115636 




0.87 


24.8000 


2.8874 


2.8735 


0.5 



2.7 Galaxy redshifts 

In star-forming galaxies such as those presented here, the 
observed Lya emission is redshifted relative to the intrinsic 
galaxy redshift, while the interstellar abs orption lines are 
blue-shifted (see e.g. IShaplev et all 120031 '). In order to im- 
prove the accuracy of the galaxy redshifts, we must account 
for these offsets. We do this us ing one of following 3 exp res- 
sions, following the method of lAdelberger et al.l (|2005aT ). 

z = z + 0.070Az - 0.0010(0-2.7) (1) 
z = z em - 0.0033 - 0.0050(z cm - 2.7) (2) 
z' = z abs - 0.0022 - 0.0015(z abs - 2.7) (3) 



where z' is the intrinsic galaxy redshift, z em is the redshift 
of the Lya emission line, z a bs is the redshift of the ISM 
absorption lines, z = (z em + z a b s )/2 and Az — Zem Z&bs . 
We note tha t we d o not use the more recent relations of 
ISteidel et all ( |201Ch : iRakic et all |201ll) as these are based 
on lower redshift BX-selected galaxies. 

Eqn. [1] relies upon having a reasonably robust measure- 
ment of both the Lya emission and the ISM features; we 
therefore use it only for those LBGs with high-quality red- 
shift IDs (Q ^ 0.6), as ISM features are only weakly detected 
in the lower-quality LBGs. For the low-Q sources we instead 
use Eqn. [2] which determines an intrinsic redshift based on 
the Lya emission line alone. 

For the LBGs showing no Lya emission the intrinsic 



redshift is calculated based on the ISM absorption features 
and using E qn. [3l with z a b a = ZL ya = zism- 

IPaper J estimated the errors on the LBG redshifts us- 
ing simulated spectra. Here, we extend the investigation into 
the redshift errors in our survey by using duplicate red- 
shift measurements. The fields presented here, particularly 
Q2348, were designed with overlapping regions and conse- 
quently there are some LBG candidates which were observed 
in more than one mask. In cases where these duplicated tar- 
gets are confirmed as LBGs, this provides two independent 
redshift measurements for the same LBG, and thus a direct 
observational test of the redshift measurement accuracy. 

Fig. ll2l shows the Az distribution for the LBGs with du- 
plicate observations, where Az = \zi — Za\ is the difference 
between the two redshift measurements. 20 objects were 
classified as LBGs in two separate observations; of these, Fig. 
[T2] indicates that 16 had fairly small errors of Az < 0.02 (of 
which 13 had very small errors of Az < 0.005), while 4 had 
considerably larger errors. In addition to these 20 objects, we 
have also searched a region o f our Q0302 field which over- 
laps with ISteidel et all (|2003l ) survey for any LBGs which 
were identified in both surveys: we find 3 such objects, and 
the redshift differences for these galaxies are also indicated 
in Fig. [H 

The standard deviation of the 20 Az values measured 
for duplicate observations in our survey is a — 0.036, corre- 
sponding to a velocity error of «2,800 km s _1 assuming a 
redshift of z — 2.8, the sample mean. However, this misrep- 
resents the true error in our redshift measurements, since it 
is skewed by the 4 sources with very high Az. These 4 values 
do not represent redshift measurement errors, rather catas- 
trophic outliers. In the cases where we find large Az val- 
ues, the error does not arise due to uncertainty in the peak 
wavelength, but in uncertainty over which spectral feature 
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Table 9. Example LBG identifications in the Q2348-011. 



ID 


R.A. 


Dec. 


u- g 


g-r 


r 


ZLyct 


ZlSM 


QlD 




(J2000) 














VLRS J235206.92-005646.70 


358.0288391 


-0.9463067 


1.30 


0.30 


24.4600 


3.1430 


2.3707 


0.5 


VLRS J235200.98-005903.ll 


358.0040894 


-0.9841969 


2.50 


0.46 


24.2400 


3.3532 


3.3435 


0.5 


VLRS J235201. 68-010002.18 


358.0069885 


-1.0006067 


1.23 


0.31 


24.9700 


2.8632 


2.8603 


0.5 


VLRS J235155.14-010104.04 


357.9797363 


-1.0177902 


1.66 


0.15 


23.8600 


2.7523 


2.7503 


0.7 


VLRS J235209.28-005535.49 


358.0386658 


-0.9265237 


1.28 


0.13 


24.0700 


2.6274 


2.6194 


0.7 


VLRS J235202.62-004747.89 


358.0109253 


-0.7966368 


2.11 


0.32 


24.2700 


3.0640 


3.0646 


0.5 



w 

CD 
CD 

D) 
CD 

33, 

< 



0.4 
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0.0 
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Figure 8. Sky distribution of LBGs and QSOs in the VLRS fields. Grey open circles show LBGs presented in Paocr % black filled circles 
show LBGs identified in the current work and cyan stars show known QSOs. 



is actually Lya. In these cases, different spectral features 
have been identified as Lya, leading to large Az. These are 
therefore better characterised as identification errors, in that 
two different solutions have been reached in the two obser- 
vations. 

For the 16 duplicated targets shown in Fig. \12b . the 
same feature has been identified as Lya and therefore the 



Az for these objects gives an indication of the measurement 
error. The standard deviation for these objects is a — 0.005, 
corresponding to Av « 380 km s _1 . 

The suggestion, therefore, is that «80% of our LBGs 
have redshift measurement errors of Aw sC 400 km s~\ while 
the other 20% may have larger errors. This problem, how- 
ever, disproportionately affects sources with an ID quality 



© 2012 RAS, MNRAS 00Q.rHf24l 



VLT LBG Redshift Survey III 13 



100 

80 

60 

40 

20 


100 

80 

60 

40 

20 


100 

80 

60 

40 

20 




□ Q>0.5 Q2359+0653 
[Q20.6 
Q>0.7 




12 3 4 



Figure 9. Redshift distributions of identified sources in each of 
the target fields. In each panel we show the distribution for the 
full sample (Q ^ 0.5), as well as for the Q ^ 0.6 (hashed) and 
Q > 0.7 (filled) subsets, where Q is the redshift quality parameter. 
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Figure 10. Redshift distribution of all identified sources in our 
5 target fields. The distribution is shown for the full sample (i.e., 
sources with Q ^ 0.5; unfilled), as well as for the Q ^ 0.6 (hashed) 
and Q ^ 0.7 (filled) subsets, where Q is the redshift quality pa- 
rameter. 



parameter Q = 0.5: of the four sources with large Az, one 
was given a quality factor of 0.5 for both redshift measure- 
ments, while the other 3 have one measurement with Q = 0.5 
and another with Q — 0.6; in the latter cases the Q = 0.6 
measurement is fairly robust while the Q = 0.5 measure- 
ment is less reliable. Therefore, the LBGs which may suffer 
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Figure 11. Redshift distribution of all identified sources our 5 
target fields, separated by initial candidate priority. Panel (a) 
shows the LBG_PRll, LBG_PRl2 and LBG_PRl3 classes, while the 
LBG_DROPs are shown separately in panel (b) for clarity. Colours 
arc as in Figs. I3H61 
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Figure 12. The distribution of redshift measurement errors, Az, 
calculated using the LBGs which were observed twice in our sur- 
vey and therefore have two independent redshift measurements. 
In panel (a) we show the full distribution, in panel (b) a close-up 
of the distribution at z < 0.025. Overplottcd as a hashed his- 
togram in panel (a) are the Az value s for 3 sources in ou r survey 
which had a redshift in the survey of lSteidel et al. (2003). 



from large errors can be excluded by removing the Q = 0.5 
LBGs from the sample. 



2.7.1 Composite spectra 

We have calculated composite spectra using the z > 2 VI- 
MOS low-resolution galaxy data. The composite spectra 
were generated by averaging over the spectra after having 
corrected the spectra for the instrument response and hav- 
ing masked skylines. In addition, each individual spectrum is 
normalised by its median flux in the rest-frame wavelength 
range 1250A ^ A ^ 1450 before being combined to form the 
composite. 
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Figure 13. Composite spectra for galaxies showing Lya in ab- 
sorption (top) and emission (bottom), with ISM absorption lines 
clearly detected. 



In Fig. [I3]we show stacked spectra for the LBGs, sepa- 
rated into those showing Lya in emission (lower panel) and 
in absorption (upper panel). These stacked spectra show 
the average ultraviolet SED of a z m 3 LBG with excel- 
lent signal-to-noise, and the quality of these spectra provide 
an indication of the robustness of our LBG identifications. 
In Fig. 1141 we show 3 separate composite spectra for sources 
classed as LBe's with quality IDs of Q = 0.5, Q = 0.6 and 
Q ^ 0.7. These spectra reflect the quality criteria set out 
in £12.51 well, with increasing quality spectra clearly showing 
increasingly high signal-to-noise in both Lya emission and 
ISM absorption features. The absorption features are only 
weakly detected in Q — 0.5 LBGs. 



2.7.2 Quasars & AG N 

We have identified 30 z > 2 AGN and QSOs in our spectro- 
scopic sample by the presence of broad emission lines. Their 
spectra are shown in Fig. [T5j All of the QSOs show broad 
Lya emission and most (at least 3 /i) also show broad C IV 
A1549, a high-ionisation emission line often used to constrain 
the dynamics and size of the broad-line r egion in QSOs as 
well as the virial mass of the black hole (Fine et alj|201(j 
and references therein). Many of the QSOs also show broad 
Ovi A1035 and Nv A1241, the latter usually being blended 
with Lya. These emission lines are indicated in each panel 
of Fig. [lU 

Several other emission lines are detected in some of the 
QSO spectra but are not marked in the figure. Ly/3 A1026 
is clearly seen in panel (d), where it may be asymmetrically 
broadened to longer wavelengths by the presence of rela- 
tively weak Ovi A1035. Panel (n) shows an emission line 
peaking at 1029A, likely suggesting a blend between Ly/3 
and O VI. 

Panels (i) and (p) both show very broad emission at 
« 1500A. These emission features cover the expected wave- 
length range of an Fein multiplet, however they consider- 
ably brighter than would be expected for these Fein lines 
(Vestergaard & Wilkes, 2001; Sigut, Pradhan & Nahar, 
2004), therefore the origin of these emission lines is unclear. 
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Figure 14. Composite spectra for galaxies classified as LBe's, 
separated by ID quality parameter Q. Those with Q = 0.5 (upper 
panel) show comparatively weak Lya emission and marginal de- 
tections of ISM absorption lines. Moving to higher quality values, 
the strength of the Lya line increases and we also see high-SNR 
absorption features. Note that the spectra are not flux calibrated. 



Ol A1304 emission appears in panels (g), (t) and (s) 
(but is difficult to see in the figure), however in all 3 cases the 
peak of the emission appears at slightly longer wavelength 
than the nominal 1304A. This may indicate a blend with 
Sill A1308, but is more likely due to the low-ionisation lines 
(which include O i) arising from a separate part of the broad 
line region to the high-io nisation lines and thus having a 
marginally higher redshift (|Constantin et al]|2002l ). Finally, 
panels (a), (c), (1), (m), (o), (r), (s) and (t) all show clear 
emission at ~ 1400A, arising from a blend of the SilV A1396 
and Ovil A1402 transitions. 



3 CLUSTERING OF z « 3 LBGS 

We now analyse the clustering of the z w 3 LBGs. As well as 
offering insights into the growth and evolution of structure in 
the Universe, we aim to measure the dynamics of the z « 3 
galaxy population, i.e. peculiar velocities and gravitational 
infall, to inform the analysis of the gas-galaxy relationship 
via LBG-Ly-a forest cross-correlation (see Paper II). 

We note that for some of the observed fields, there is 
some variation in the magnitude cuts applied in selecting 
the candidates for spectroscopic observation. The analysis 
of the clustering is therefore limited to those galaxies with 
R < 25. Aside from this magnitude cut, all galaxies with 
Q ^ 0.5 are included throughout this analysis. 

In the analysis that follows, we measure the galaxy clus- 
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Figure 15. Rost-framc VIMOS spectra for the QSOs detected in 
the 5 fields presented. The j/-axis scales differ from panel to panel, 
but the dotted line marks zero flux in each case. Dashed red lines 
indicate the wavelength of O VI, Lya, Ne V and C IV emission. 
Gaps in the spectra indicate that an artefact has been masked 
out. 
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Figure 16. Pair counts (DDs) as a function of pair separation in 
the transverse, a, and 3D Hubble, s, directions. The Keck data 
(dashed orange curves) provide better sampling at small scales, 
whilst the VLRS data (solid black curves) provides sampling at 
larger scales that is not provided by previous data. 



tical area. The redshift distribution of the random catalogue 
is set using a polynomial fit to the data. 

We focus on fitting the auto-correlation function in the 
semi-projected, w p (a), and full 2-D, £(<r, tt), forms, where 
a and tt are the separation of two galaxies transverse and 
parallel to the line-of-sight. But we shall also study the LBG 
z-space correlation function, £(s), where the signal can be 
significantly higher at large scales. 

For £(cr, tt) in particular, we also consider the com- 
bi ned sample o f the V LRS data with the Keck LBG data 
of ISteidel et~aH (|2003T ). The Keck data consists of 813 z « 3 
galaxies and offers higher sampling rates than the VLRS, but 
across smaller field sizes (ss 8'). This is illustrated in Fig. 1161 
where the solid black line shows the VLRS pair counts (DD) 
as a function of separation in the transverse direction, a 
(left hand panel), and the 3-D separation, s = y/a 2 + tt 2 
(right hand panel). In both panels the Keck pair counts are 
shown by the dashed orange line. Fig [16] shows that the 
VLRS pair counts in both the transverse and 3-D distance 
are significantly higher than in the previous Keck sample at 
a,s > 10h _1 Mpc. 



3.1 Semi-projected correlation function, w p (a) 

We first estimate the LBG clustering using the semi- 
projected correlation function, w p (a). This gives the clus- 
tering at fixed transverse separation, a, integrated over line- 
of-sight distance, tt, approximately independent of the effect 
of peculiar velocities and is given by: 



tering as a function of galaxy-galaxy separation using the 
Landy-Szalay estimator: 

= (DD(x))~2(DR(x)) + (RR(x)) (4) 
(RR(x)) 

where £(x) is the clustering as a function of separation x, 
DD(x) is the number of galaxy-galaxy pairs at that sepa- 
ration, DR(i) is the number of galaxy-random pairs and 
RR(x) is the number of random-random pairs. This is es- 
timated using a random catalogue that consists of 20 x as 
many random points as data points and that covers an iden- 



w p (a) = 2 Z(a,TT)dTT (5) 
Jo 

where £(cr, tt) is the 2-D auto-correlation function. We inte- 
grate £(cr, tt) over the range < tt < 7r max (cr), where 7T ma x(c) 
is given by the maximum of 1000(1 + z)/ H(z) and 15cr at 
a given sky separation a (co nsistent with lAdelberger et all 
l2003l : Ida Angela et al.ll2005al ) . 

In the calculation of w p (a), we make a c orrection for 
the effect of 'slit collisions', following IPaper A Any object 
observed with VIMOS takes up an area on the detector of 
at least 40 x 640 pixels ( jjgXTJ , corresponding to 8" x 130" 
on-sky. Other candidates lying within this area can there- 
fore not be observed (unless the area is revisited), and as a 
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result, pairs of LBGs at small separations are systematically 
missed by our survey. This effect will redu ce the m easured 
LBG auto-correlation at small separations. iPaper J quanti- 
fied this effect by comparing the angular auto-correlation 
function of photometrically selected LBG candidates and 
spectroscopically observed candidates. Using their result, we 
correct for this effect in our LBG survey by weighting DD 
pairs at 6 < 8' according to the weighting factor given by 



W(6) 



1 



1 - (0.0738 x 



(6) 



where 9 is the angular separation in arcminutes. 

In addition to the slit collision correction, a further cor- 
rection - the integral constraint - is required to compensate 
for the effect of the limited field sizes. For this we follow 
the commonly used approach of using the random-random 
pair distributions, which have been constructed to match the 
survey geometry, to determine the magnitude of the integral 



ber of authors (e.g. Groth & Peebles! 


1977; Peebles 


198d; 


IRoche et all 19931; 


Baugh et al.l 1 19961: 


Phleps et al.l 


20061). 


with Phleps et al.l ( 


2006T) in particular providing: a detailed 



discussion in relation to the projected correlation function, 
and we provide a brief description of the calculation here. 

The measured correlation function is given by the true 
correlation function minus the integral constraint, X: 



Wm(a-) = Wt(a) — X 



(7) 



Assuming a power-law form for the the real-space clus- 
tering, the true projected clustering is fit by: 



w p (a) = CVqU 



7 1—7 



(8) 



where ro is the real-space clustering length and 7 is the 
slope of the real-space clustering function, £(r), which is 
characterised by a power-law of the form: 



The factor C is given by: 



C = 



r(i)iW) 



(9) 



(10) 



where F is the Gamma function. Given this framework, the 
integral constraint can be estimated from the mean of the 
random-random pair counts, (RR) , and the slope of the cor- 
relation function, such that: 



X 



T,(RR(a))a 
E (RR(a)} 



1-7 



(11) 



Quantifying errors on the auto-correlation function has 
been performed using Poissonian, jack-knife and random re- 
alisation error estimates. The Poisson errors are given by: 



(i + O 

V(DD) /2 



(12) 
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Figure 17. Ratio of the random mock (red dashed histogram) 
and the jack-knife (blue dot-dashed histogram) error estimates to 
the Poisson errors as a function of separation for the projected 
correlation function (left panel) and the redshift space autocorre- 
lation function (right panel). The results are shown for the VLRS 
sample only. 



The jack-knife errors were computed by splitting the 
data into individual fields, with the large fields (i.e. HE0940 
and Q2348) being split into two fields each. We therefore 
have 11 different jack-knife realisations with a single field 
(or half-field) being excluded in each realisation. 

The random realisation error estimates incorporate 100 
random catalogues with the same number of objects as the 
real data. We then calculate the correlation function using 
these random realisations to calculate the {DD) pairs and 
take the standard deviation of the results as the uncertainty 
on the measurement. 

In Fig. 1171 we compare the above error estimates for 
w p (a) and the redshift space clustering function, £(s) (see 
section [3^2]) . showing the ratio of the jack-knife and random 
realisation methods to the Poisson result. The estimates are 
consistent with each other over scales from « 1 — 25 /i" 1 Mpc. 
In what follows, we therefore use the Poisson estimates at 
separations of < 25 /i -1 Mpc and jack-knife estimates at sep- 
arations > 25 ft. _1 Mpc. 

The projected auto-correlation function for the VLRS 
sample (black cir cles), the Keck sample (orange squares; 
ISteidel et al.ll2003T ) and the two combined (cyan triangles) is 
shown in Fig [18] The plotted data include the integral con- 
straint correction, which we estimated as I wp — 5.29 and 
Imp = 7.18 for the VLRS and Keck data respectively. 

Based on this we estimate a clustering length of the 
entire VLRS sample of ro = 3.32±0.41 /i _1 Mpc (comoving) 
with a slope of 7 = 1.53 ± 0.13. The Keck result on it's 
own gives a result of ro = 3.98 ± 0.32 h~ 1 Mpc with 7 = 
1.58 ± 0.13, whilst the combined VLRS+Keck data gives 
r = 3.75 ± 0.24 h^Mpc with a slope of 7 = 1.55 ± 0.09. 
These ro results are comparable to the clustering of star- 
forming galaxies at lower redshifts (e.g. iBlake et al.ll2009l : 
iBielbv et af]|2010D . 

Comparing t o other measurements o f the z ~ 3 LBG 
clustering length, iGiavalisco fc Dickinson! (|200ff ) measured 
r = 5.0 ± 0.7 /i _1 ]VIpc for Rab < 25 LBGs, but for a 
relatively s mall number of ga l axies (« 400). Building on 
that sample. I Adelberger et all (|2003h measured ro = 3.96 ± 
0.29 h _1 Mpc with aslope of 7 = 1.55±0.15 at Rab < 25.5. 
We note that with the sa me sample, but a different method, 
lAdelberger et all (|2005br i found a highe r clustering strength 
of r = 4.5±0.6 /i _1 ]Vipc. Subsequently. ICooke etafl (|2006t ) 
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Figure 18. Projected auto-correlation function, Wp(a)/cr, for the 
VLT LBG sample (black circles). The result for the 2.5 < z < 3.5 
Keck LBG sample is shown by the orange squares and error bars 
for both results are estimated using Poisson errors. The best fit 
power laws for both the VLT (solid black line) and Keck (dashed 
orange line) are shown, with clustering amplitudes of ro = 3.32 ± 
0.41 fc-iMpc and r = 3.98 ± 0.32 h^Mpc, with slopes of 7 = 
1.53 ± 0.13 and 7 = 1.58 ± 0.13 respectively. The cyan triangles 
show the combined result and the dotted cyan line the best fit to 
this of r = 3.75 ± 0.24 h^Mpc with a slope of 7 = 1.55 ± 0.09. 




Figure 19. The two-dimensional auto-correlation function, 
£(<T, tt) results for the VLRS (left) and Keck (right) data sam- 
ples individually. The shaded contour map and dashed contour 
lines show the measured §(<r, tt) in each case. 



3.2 2D Auto-Correlation Function, £(<r, rr) 

As discussed above, integrating along the redshift/line-of- 
sight direction leaves the semi-projected correlation func- 
tion, w p (a), independent of the effects of galaxy peculiar 
motions. We now attempt to fit the full 2D correlation func- 
tion, £(cr, 7r), to retrieve the kinematics of the galaxy popu- 
lation and to make new estimates of ro. 

As before, we use the Landy-Szalay estimator to calcu- 
late the correlation function but now as a function of both 
transverse separation, a, and line-of-sight separation, n. We 
use the same random catalogues matching the survey fields 
as used for the calculation of the projected correlation func- 
tion. We again calculate the integral constraint for the data 
sets using the random catalogues via: 



measured the clustering of z ~ 3 LBGs in fields around 
damped Lya absorbers and found a lower clustering strength 
of r = 2.65±0.48 fr _1 Mpc with a slope of 7= 1.55±0.40 at 
Rab < 25. Overall, our result therefore appears consistent 
wit h other measurements, a lthou g h marginally lowe r than 
the Giavalisco fc Dickinson! |200ll ); lAdelberger et all |2003l . 
l2005bl ) results, which are all based on the same - or a subset 
of the same - sample. 

The above estimates are based on spectroscopically 
confirmed samples and a number of clustering measure- 
ments exist based o n photometric samples. For example 
iFoucaud et~ai1 l|2003l ) measur ed r = 5.9±0.5 /i _1 Mpc from 
the angular correlation function of Rab < 24.5 LBGs from 
the C anada-France Deep Fields Survey ()McCracken et all 
l200ll ). a higher ro than the spectroscopic samples, but 
also a significantly b r ighter magnitude cut. Addition- 
ally, lAdelberger et all (|2005bl ) measured w(6) for pho- 
tometrically selected LBGs and found ro = 4.0 ± 
0.6 fc~ 1 Mpc for R < 25.5 LBGs, consistent with our results. 
IHildebrandt et all (|2007h measured r = 4.8 ± 0.3 h~ 1 Mpc 
for 22.5 < -Rvega < 25.5 galax i es in the GaBoDS data. Sub- 
sequently, IHildebrandt et all (|2009l ) measured ro = 4.25 ± 
0.13 fc _1 Mpc for CFH TLS LBGs at tab < 2 5 and using 
photo-z from HYPERZ (|Bolzonella et al.ll200ol ). Again, our 
estimate for ro is marginally lower than these measurements, 
however as noted the photometric selections will be prone 
to a certain level of uncertainty caused by contamination 
corrections. 



1 _ £(AR(s))s^ 



£ {RR(s)) 



(13) 



where s — \f[ry 2 + tt 2 ). This gives values of = 0.036 and 
1% = 0.064 for the VLRS and Keck data samples respec- 
tively. 

Fig. [19] shows the result for the VLRS data (left-hand 
panel), which provides a greater handle on the large scale 
(s > 10ft _1 Mpc) clustering, the Keck data (right panel), 
which provides greater sampling on small scales. Fig. [20] 
shows the VLRS and Keck results combined. In each case, 
f(<7, 7r) was calculated in linear 2 /i~ 1 Mpc bins and subse- 
quently smoothed with a fwhm of 2 fe _1 Mpc. 

For both the VLRS and Keck samples, we see the 'finger 
of god' effect at small a scales in which the clustering power 
is extended in the tv direction. This effect is a combination 
of galaxy peculiar velocities and measurement errors on the 
galaxy redshifts. In addition, in the VLRS a flattening of 
the clustering measurement at large scales is evident, which 
is caused by dynamical infall of galaxies. 

We now fit models of the clustering to these results, 
initially assuming a single power-law for £(r) and allowing 
ro and the kinematical parameters to vary. We take the ro 
and 7 estimates from the w p (a) fit as the starting point in 
fitting the 2D clustering. The kinematics are characterised 
by two parameters: the velocity dispersion in the line of sight 
direction \J (tof) and the infall parameter, /3. The model 
we use incorporating the gal axy k inematic s is described in 
full by lHawkins et"al] (|2003l ) and IPaper J . For the VLRS 
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Figure 20. The two-dimensional auto-correlation function, £(cr, 7r) results for the VLRS sample (left) and the signal-to- noise on the 
result (right). The shaded contour map and dashed contour lines show the measured £(<r, tt), whilst the solid contour lines give the best 
fitting model. 



and the combined samples we fit over the range 1.0 ^ s ^ 
25 /i -1 Mpc, whilst for the Keck data by itself we limit the 
fit to the scales 1.0 ^ s ^ 15 /t _1 Mpc (note that the largest 
single field available in the Keck data is « 15 ft - Mpc). 

For the two samples individually, we find that it is 
difficult to place reasonable constraints on both the ve- 
locity dispersion and the infall together. With the VLRS 
data (over the range 1 ^ i< J ,' lr ) ^ 25 ft, - Mpc), we find 
/3(z = 3) = O.Sto.l and y/twf) = lTOOi^kms" 1 , the 
low signal-to-noise on small scales limiting the fit accu- 
racy. We experimented with adding a uniform error distri- 
bution out to ±12000kms _1 to the Gaussian velocity dis- 
persion (c.f. Fig. I12p but this made little difference in the 
a, tv range fitted. Fitting the Keck data gives best fit values 
of /3(z = 3) = 0.85±g |g and y^w f) = 700 ± 220km S - 1 . We 
note that Ida Angela et al.l (|2005al ) performed a similar fit to 
the Keck data for f3(z = 3), but kept a constant velocity 
dispersion of \J (tof) = 400kms -1 , finding a value for the 
infall parameter of j3(z = 3) = 0.15!q'j5. By also setting 
the velocity dispersion to a value of 400kms~ 1 , we find tha t 
we retrieve a comparable result to Ida Angela et all (|2005al ). 
highlighting the degeneracy between the velocity dispersion 
and the infall parameter. 

Ultimately, fitting the VLRS £(cr, n) is hindered by a 
lack of signal-to-noise on small scales, whilst the fit to the 
Keck data is hindered by the small size of the fields. We thus 
combine the two datasets and fit the full LBG sample in the 
same manner as with the individual samples. The fit is per- 
formed in the range 1 < s < 25 /i _1 Mpc and we allow the ve- 
locity dispersion and the infall parameter to vary. The result- 
ing fit gives a velocity dispersion of ^JTwT) = 540t2ookms -1 
and an infall parameter of /3(z = 3) = 0.32 ± 0.20. We show 
the x 2 contours for the fit in the /3(z = 3) — -J (w'i) plane 
in Fig. [21] (the contours represent the la, 2a and 3a confi- 
dence limits). From this figure, the degeneracy can be seen 
between \J and j3, where increasing ft similarly increases 
the best fit velocity dispersion. The best fitting results are 
plotted over the contour maps of the £(cr, 7r) measurements 
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Figure 21. Fitting contours for peculiar velocity and bulk in- 
flow based on the combined VLRS+Keck £((X, tt). The best fit- 
ting result is given by /3(z = 3) = 0.32 ± 0.20 and y/{w'£) = 
540l^kms- 1 . 

in Fig. [501 (dashed contours) . As with the data, we see the 
finger-of-god and large scale flattening effects in the fitted 
models. 

3.3 Redshift-space correlation function, £(s). 

In order to check the consistency of our measurements, 
we now compare the model fit obtained from w p (a) and 
£ (a, tt) to the measured redshift-space auto-correlation func- 
tion £(s). Again we use the Landy-Szalay estimator and 
quote errors based on Poisson estimates. The £(s) results 
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Figure 22. Redshift space correlation function, for the 

VLRS (open blue circles), Keck (orange squares) and combined 
(black triangles) samples. The short and long dashed curves show 
single power law model fits to the combined sample based on fits 
to uip(o-) and f;(cr, 7r). The dash-dot curve gives a double power 
law model fit, whilst the dotted curve at large scales shows the 
predicted ACDM clustering. 



for the VLRS, Keck and combined LBG samples are shown 
in Fig. I22I We also plot the single power law estimate of the 
intrinsic clustering from our fits to the VLT+Keck £(cr, n) 
(dotted line) and the result of this power-law after applying 
the best fit values for /3 and \J (uif ) (dashed line). The fi- 
nal fit incorporating the galaxy dynamics is marginally low 
compared to the data-points, but is easily consistent within 
the error bars. 

Our measurements of /3 and v(wi) are consistent with 
the prev ious measureme nts usin g the first VLT dataset 
|Paper J ). As discussed in lPaperll . the median measurement 
error on the galaxy lines on the VLT VIMOS spectra is 
« SGOkms" 1 . In addition, an uncertainty of ~ 200kms 
is introduced by the transformation from outflow redshifts 
to intrinsic galaxy redshifts using Eqns [TJ [2] and [3] The 
final contribution to the velocity dispersion is from the 
intrinsic pec uliar velocities of t he galaxies. Using the GIMIC 
simulations (|Crain et alj|2009l ) we have analysed the mean 
velocity dispersion of LBG-like galaxies and find a value of 
~ lOOkms -1 . Combining these three elements in quadrature, 
we would expect a pairwise velocity dispersion of y/ ~ 
y/2 x \Z(360kms -1 ) 2 + (200kms" 1 ) 2 + (lOOkms" 1 ) 2 « 
600kms _1 . This is within the la error contours for a 
slightly higher value of /3 in Fig. 1211 This value is also 
reasonably consistent with the VLRS £(s) estimate (see 
Fig. [22]). 

A ACDM £(r) with standard values for the cosmological 
parameters (flat space, tt m = 0.27, H = 70kms~ 1 Mpc _1 
and flbh 2 = 0.022) is also compared to the combined data 
in Fig. [55] where the linear model has been normalised to 
the combined £(s) at 10 h~ Mpc. There are currently some 
claims that non-Gaussianity is d etected at z ~ 1 in NRAO 
VLA Sky Survey radio source dXia et alj 120101) and Lu- 
minous Red Galaxy (LRG) datasets ([Thomas et alj l201ll ; 
ISawangwit et al.ll201ll ; lNikoloudakis et alJIln PretJ) . The ev- 



idence generally comes comes from detecting large scale ex- 
cess power via flatter slopes for angular correlation func- 
tions. Since non-Gaussianity is easier to detect at high red- 
shift this motivates checking the LBG £(s) for an excess. 
We have already noted that the slope from w p and £(e>", n) 
at 7 = 1.55 is much flatter than the canonical 7 = 1.8. 
This slope is also flatter than the z ~ 1 L RG large-scale 
w(9) slope of iNikoloudakis et"al] (|ln Prepl) . We see that 
the VLRS does give reasonably accurate measurements for 
25 < s < 100 fe _1 Mpc and that the observed LBG £(s) 
shows a ~ 2.5<r excess over the ACDM model in this range. 
Even when the marginally smaller integral constraint for the 
ACDM model is assumed the discrepancy remains at m 2a. 
We conclude that there is some evidence for an excess over 
the standard ACDM model but independent LBG data is 
needed to confirm this on the basis of the redshift space 
corre lation fu nction. The statistical error on the LBG w(9) 
from iPaper il is smaller but the flat power-law here is only 
seen to 9 — 10' or r — 13— 14 /i _1 Mpc and this is not enough 
to decide the issue. 



3.4 Double power-law correlation function models 

We next look to see if a more complicated model than a 
powe r-law for £(r) is required. This is motivated firstly be- 
cause IPaper J noted that there was an increase in the slope 
at ~ 1 /i _1 Mpc in the LBG angular auto-correlation func- 
tion, w(9), suggestive of the split between the 1-halo and 
2-halo terms in the halo model of clustering. Although this 
result is uncertain due to quite significant low redshift con- 
tamination corrections, such features have been seen in lower 
reds hift galaxy samples, particularly for LR Gs at z ~ 0.5 
(e.g. lRoss et al.; 2007: Sa wangwit et al.ll201ll ). Given the im- 
proved power of the VLRS, it is interesting to see if there is 
any evidence of a change in the slope at small scales in £(s) 
and w p (a) in our z ~ 3 LBG sample. 

We therefore show in Fig. [18] a double power-la w model 
for Wp with the same power-law slopes as fitted by IPapeFll 
to the LBG w(9). We have reduced the amplitude by ~ 
20% to match approximately the large-scale amplitude fitted 
to the VLRS and Keck combined data. This is within the 
systematic uncertainties of the w(9) measurement. Although 
certainly not required by the w p data this double power- 
law cannot be rejected by the combined w p data, giving a 
reduced x 2 of 1.77 (marginally smaller than the reduced x 2 
obtained for a single power law of 1.84). 

In Fig. [22] we now compare to £(s) the same double 
power-law w{9) model with the « 20% reduced amplitude. 
Again with a velocity dispersion of 540kms~ 1 and /3 = 0.32 
we see that the model cannot be rejected by the data. We 
note that if we use a (DD) / (DR) estimator the VLRS £(s) 
result shows increased power at large scales and the flatter 
slope of the double power-law model here provides a better 
fit. 

We note that other authors have also reported a turn- 
up in the clustering at small scales in h igh redshift galaxy 
samples. For instance lOuchi et al.l (|2005l ) reports that z = 4 
LBG w(9) shows a steepening below « 0.2 /i _1 Mpc or « 10" 
at z = 4. If both results are unaffected by contamination 
then it could argue for an evolutionary growth in this break 
scale between z — 4 and z — 3. 

Certainly there is plenty to motivate expanding surveys 
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Figure 23. Fitting contours for the mass density and bulk in- 
flow based on the combined VLRS+Kcck £(ct, 7t). The shaded 
region gives the result from the VLT+Keck data sample (as- 
suming y/ (w%) = 540kms~ 1 . ro = 3.75 h~ 1 Mpc and 7 = 1.55) 
with the dotted contour lines giving the 1, 2 and 3a uncertain- 
ties. The results from this fit are /3(z = 3) = 0.38^q and 
n m (z = 0) = 0.08ioo8- The dashed red lines show the 1, 2 
and 3cr constraints given by evolving the 2dFGRS measurements 
as described in the text. The solid black contours give the com- 
bination of the two and give a result of P(z = 3) = 0.38±o;Jg and 
n m (z = 0) = 0.30±g;?|. 



of 0. m (z = 0) = 0.081q'o8 and on the infall parameter of 
0{z = 3) = 0.38i°;i». 

As discussed bv lda Angela et al] (|2005al lbh and lPaper J . 
by combining these results with a measurement of the am- 
plitude of the mass clustering at low redshift, it is possible 
to break the degeneracy between Q, m and (3. To do this , 
we use the 2dFGRS measurements of lHawkins et al' I (|2003h . 
who found r = 5.0 /i _1 Mpc, 7 = 1.8 and /3(z = 0.11) = 
0.49 ± 0.09. We then determine the z — 0.11 underlying 
dark matter clustering amplitude from these parameter con- 
straints and evolve this to z = 3 for test cosmology range 
of Q m (z = 0) =0—1. Finally, we estimate the bias for 
the LBG population from this extrapolated matter cluster- 
ing amplitude and the measured LBG clustering properties 
we have presented in earlier sections. The constraints on /3 
using this method over a range of assumed fl m values are 
given by the red dashed contours in Fig 1231 By combining 
these with the original constraints from £(cr, 7r), we find a 
result of P(z = 3) = 0.38±°;^ and Q m {z = 0) = 0.30±J^g. 

Across these analyses, we have consistently found a 
value for the infall parameter of f3(z — 3) m 0.36 — 0.40. 
flm is somewhat less well constrained, but remains consis- 
tent with ACDM. The measurements of P(z = 3) presented 
here are consistent with our previous measurement from 
iPaper J of /3(z = 3) = 0.48 ± .17, whilst being somew hat 
higher than the result found bv lda Angela et al.l (|2005al ) of 
/3(z = 3) = 0.15^0 is- We note that the latter assumes a 
fixed velocity dispersion of J {wV} = 400kms~ 1 and given 
the degeneracy between /3 and y{wf) their lower estimate 
of P is not surprising. 



to make more accurate measurements of both the angular 
and redshift survey correlation functions at these redshifts. 
Below the break scale is of extreme interest for single halo 
galaxy formation models and at large scales the interest is 
in looking for a flattening of the correlation function slope 
due to the presence of primordial non-Gaussianity. 



3.5 Estimating fi m and the growth rate 

3.5.1 The mass density of the Universe 

We now look at the cosmological results afforded by the 
z « 3 LBG clust ering and dynam ics. As discussed by 
iHovle et al] (j2002l ): Ida Angela et al] (|2005bl ). it is, in prin- 
ciple, possible to constrain the matter density Q m (z = 0) 
from the measurement of £(<r, n). Effectively, the elongation 
of £(cr, 7r) along the line of sight increases with increasing 
values of Q. m {z = 0). However, increased values of /3 lead to 
a flattening of £(<r, ir) along the line of sight. These effects 
combined lead to a degeneracy in determining Q, m from the 
galaxy clustering alone. 

In previous sections, we have studied the galaxy dynam- 
ics assuming a cosmology with Q m (z = 0) = 0.3. We now 
fit the £(cr, 7r) result with this assumed cosmology, but now 
with a constant peculiar velocity of \J (wf ) = 540kms _1 and 
fitting for Q m (z = 0) and /3. The result is shown by the la, 
2a and 3a contours (solid) in Fig. 1231 Based on just the 
z = 3 galaxy clustering, we find results for the mass density 



3.5.2 Growth rate results compared 

In a similar fashion to the previous section, we may use 
the clustering and dynamics measurements to constrain the 
growth of structure via the growth parameter, f(z). The 
bulk motion and the clustering bias, b, of a galaxy popula- 
tion are related to the growth parameter via: 

p= m K w 55 (14) 

Therefore, given measurements of the bias and dynam- 
ical infall of a range of galaxy populations across cosmic 
time, constr a ints o n the growth of structure can be made. 
iGuzzo et al.l (|2008l ) presented the results of such an analysis 
based on the VLT VIMOS Deep Survey (VVDS), showing 
values for /(z) extracted from a number of galaxy surveys 
up to a redshift of 2 w 0.8. Here we add the z ~ 3 result 
from our survey. We present measurements in terms of both 
f(z) and fas, where fas is intended to give a measurement 
which is less dependent on the c osmology assumed for the 
calculation of the clustering (e.g. ISong fc Percivalll200^ ). 

We have already calculated the infall parameter and 
take the value (/3 = 0.32) obtained via fitting the velocity 
dispersion and /3 in a ACDM cosmology with Q m — 0.3 and 
Q A = 0.7 (Fig. mi). 

The bias can be calculated directly from the clustering 
measurements by using the volume averaged clustering: 
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Figure 24. The evolution of the growth factor based on available 
z < 1 survey observations and the constraint from the VLRS data 
(solid blue circle). 



Us) 



(15) 



where f s (8) is the volume averaged correlation function at 

s < 8 h~ 1 Mpc for the galaxy population and £ DM (8) is the 
same, but for the underlying dark matter distribution. The 
volume averaged clustering is calculated from the clustering 
using: 



£(r)dr 



(16) 



In addition, a measure of the dark matter clustering is 
required in order to estimate of the bias of the galaxy popu- 
lation and we calculate this using the CAMB so ftware incor- 
porat ing the HALOFIT model of non-linearities (|Smith et al.l 
l2003h . 

Using the previously determined best fit parameters of 
ro = 3.75 /i -1 Mpc and 7 = 1.55, we evaluate the galaxy bias 
based on a single power-law, finding a bias for the LBGs 
of b = 2.59 ± 0.11. Combining this with our measurement 
of the infall parameter of j3 — 0.32 gives a value for the 
growth parameter based on the combined LBG sample of 
f(z = 3) = 0.83 ±0.46. 

We present the f(z = 3) result (filled blue circle) in 
the top panel of Fig. [24] alongside a number of other low- 
redshift measurements. In ord er of ascending redsh ift, the 
star shows the measurement of lTurnbull et all (|2012l ) based 
on local supernovae measurements, the dia mond shows the 
result based on the 2dFGRS presented by I Hawkins et al.l 



|2003l ). the red triangle shows t he SDSS result based on 
LRGs from lTegmark et all (|2006t ). the red square shows the 
2SLAQ result als o estimated from the LRG population of 
iRoss et al l J2007I) and the black circle shows the VVDS re- 
sult of Guzzo et alj |2008l ). For completeness these are also 
summarised in Table [T2l 

In the middle panel of Fig. 1241 we also plot the evolution 
of f(z) based on the assumed ACDM cosmology (dashed 
line), where f(z) = Q m (z) ' 55 . The low redshift data points 
are all consistent with the assumed cosmology at the ~ la 
level and at z = 3, the model cosmology is again consistent 
with the data. We note again that the observations them- 
selves depend on the assumed cosmology via crg(z) and so 
to test the Q m = 1 cosmology we adjust the observed values 
of f(z) for the effects of different cosmology in eq. 1151 ac- 
cording to the methods set out by Ida Angela et all (|2005al ). 
Also assuming that f3 is approximately independent of the 
assumed cosmology, we see that the fi m = 1 z-independent 
growth rate is apparently rejected by the data. However, if 
the bias is allowed to float rather than just fit the lowest 
redshift point then the model may only be rejected at the 
1 — 2(7 level, consistent with the conclusions from Fig. 1231 

If we now consider /erg, the observations are now inde- 
pendent of the assumed cosmology, at least given again the 
assumption that the observed j3 is approximately cosmol- 
ogy indpendent. Each of the observational measurements is 
again plotted in the top panel of Fig. 1241 but now in terms 
of /erg. We now plot three test cosmologies for compari- 
son, the ACDM used in the top panel (dashed line), plus 
an Einstein-de-Sitter model (fl m = 1, Q,a = 0, dotted line) 
and an open Universe without a cosmological constant and 
a mass density of Q m — 0.3 (dot-dash line). For each model 
we incorporate a factor (c a t) to correct for the cosmology 
assumed in the measurement of the clustering observations 
being different from the test cosmology. Each model is thus 
given by: 



(/o"g)a = P(Tg,a = 



/Cat 



(17) 

where an index of t denotes a parameter calculated in the 
test cosmology and an index of a denotes a parameter cal- 
culated in the assumed cosmology (i.e. ACDM). We nor- 
malise £7g to 0.8 at z = and o g is effectively erg measured 
for the galaxy population. If we assume a power-law form 
for the clustering with a slope of 7 = 1.8, th en following 
iBallinger et all (|l996l ): Ida" Angela et all (|2005al ) c at is given 
by: 



Cat = 



Bt\ M_ 

Ba) A a 



2/3 



(18) 



with A and B as defined by Ida Angela et al.l (|2005af ). We 

are assuming here that j3 is independent of cosmology, a 
reasonable approximation when il m > 0.1 

With the models corrected to account for differences 
between the cosmology assumed for the observations and 
the model cosmologies, the observations now provide clearer 
tests on the models. We note that the observations (except- 
ing the SDSS point) assume the Q m = 0.3/^a = 0.7 cosmol- 
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Table 12. Summary of growth parameter results from the literature. 



Survey 


z 


b 


P 


/ 


/C8 


First Amendment SNe 1 


0.025 








0.40 ± 0.07 


2dFGRS 2 


0.11 


1.15 ± 0.06 


0.47 ±0.18 


0.54 ±0.21 


0.41 ±0.16 


SDSS LRGs 3 


0.35 


2.13 4 


0.31 ±0.04 


0.66 ±0.07 


0.44 ± 0.05 


2SLAQ LRGs 5 


0.55 


2.02 ± 0.10 


0.45 ±0.05 


0.91 ±0.11 


0.55 ±0.07 


VVDS 6 


0.77 


1.30 ±0.10 


0.70 ±0.26 


0.91 ±0.07 


0.50 ±0.04 


VLRS+Keck 


2.90 


2.59 ±0.11 


0.32 ±0.20 


0.83 ± 0.46 


0.22 ±0.12 



1 | Turnbull et alj <2012f> 2 lHawkins et alj [120031 1 3 lTegmark et alj feOOrj l 
4 Song fc Percivall j200sh 5 lRoss et alj ||2007T| 6 lGuzzo et alj J2008h 



ogy and so there is no change in the relationship between 
the observations and the ACDM model between the top and 
bottom panels. Thus the VLRS datapoint shows the same 
level of consistency with the ACDM model for fas and f(z). 

At redshifts of z < 1, we see that flat and open cos- 
mologies (i.e. with and without a cosmological constant) are 
poorly distinguished by the available observations. At z = 3, 
we find that the VLRS data rejects the open cosmology with 
fim = 0.3 at the « 2a level. The Einstein-de Sitter cosmol- 
ogy is apparently rejected by combining the z < 1 and z « 3 
observations. But again if the normalisation of the model is 
allowed to float rather than be fixed on the low-redshift SNe 
observation, the model still fits the data with a reduced \ 2 
of 2.7. 



4 CONCLUSIONS 

We have presented the widest area spectroscopic survey of 
galaxies thus far in the redshift range 2 < z < 3.5, based on 
observations with the VLT VIMOS instrument . This paper 
adds to the initial dataset of lBielbv et alj (|201ll ) , where data 
in five 0.5° x 0.5° fields were presented. Here we add a 
further four new fields, each with deep optical imaging over 
an area of w 0.5° x 0.5° in three cases and a full 1° x 1° 
in the fourth field. I n additi on, we have extended one of 
the original fields of iPaper J to 1° x 1° from the original 
0.5° x 0.5°. In total therefore, we now have « 4 deg 2 of 
optical imaging with a minimum of three bands in each field 
incorporating U, B and R or equivalents. 

In total, the survey now consists of 2,148 z > 2 galaxies. 
The properties of the full sample have been presented here 
with redshift and magnitude distributions as well as exam- 
ple and composite spectra. The mean redshift of our z > 2 
galaxy dataset is z — 2.78. In addition, we detect 30 AGN 
or quasars, w 800 low-redshift galaxies and «130 Galactic 
stars. Using the z > 2 galaxy dataset, we have conducted an 
analysis of the galaxy clustering at z ~ 3. Using the semi- 
projected correlation function, we have measured a galaxy 
clustering length of ro = 3.32 ± 0.41 /i~ Mpc with a slope 
of 7 = 1.53 ± 0.26, assuming a power-law form to £(r). We 
have also combined the VLRS sample with the Keck LBG 
sample of ISteidel et al.1 |2003|). which provides greater sta- 
tistical power on small scales (i.e. s < 2 ft _1 Mpc) than the 
VLRS but does not provide the coverage of the VLRS at 
larger scales (i.e. s > 8 /i _1 Mpc). For the combined sample 
we measure a clustering length of ro = 3.75 ± 0.24 /i~ Mpc, 
with a slope of 7 = 1.55 ± 0.09. 



We have shown that the LBG correlation functions con- 
sistently show slopes that are significantly flatter than the 
canonical 7 = 1.8 observed at low redshift. Indeed, the mea- 
sured slopes of 7 = 1.55 are flatter than in some z fa 1 galaxy 
and radio-source correlation functions that have been inter- 
preted as showing evidence for primordial non-Gaussianity 
JXia et al.ll20ld;lThomas et al.ll20l"ll ; ISawangwit et alj|201ll ; 
iNikoloudakis etal Jlln Prenl) . Non-Gaussianity is expected to 
be easier to detect at large-scales and high redshift. We have 
therefore checked whether a standard ACDM model is con- 
sistent with the form of the VLRS f (s) in particular in the 
regime 10 < s < 50 7i~ Mpc. We found that there is ev- 
idence that the LBGs are showing more large-scale power 
than the standard model in this regime but only at « 2cr. 
More studies of LBG clustering at large-scales are clearly 
needed to check these results. 

In addition to the ID clustering analyses, we have also 
investigated the 2D correlation function and the imprints of 
galaxy dynamics on the clustering. We find that the 2D clus- 
tering for the VLRS+Keck LBG sample is well fit by a model 
based on a power-law fit with a clustering length of ro = 
3.32 h~ Mpc, a large scale infall parameter of /3 = 0.32±0.20 
and a velocity dispersion of J = 540^2 O ! o ) kms~ 1 , over 
a range of 1 < s < 25 /i -1 Mpc. We have shown that this 
result is consistent with the model for the redshift space cor- 
relation function, £(s), measured for the combined sample. 

Finally, we use the 2D galaxy clustering results to deter- 
mine the matter density parameter and the growth param- 
eter. Using the previously constrained form for the cluster- 
ing and galaxy velocity dispersion, we fit the 2D correlation 
function for the matter density, Q m - We find an acceptable 



range in the matter density of Q, m {z = 0) = 0.08 
an infall parameter of /3(z = 3) = 0.38 



+0.15\ 



+0.22 
-0.08 



(with 
We add a fur- 



ther constraint provided by the 2dFGRS low-redshift clus- 
tering measurements, which gives 0, m (z = 0) = 0.30!to'i8 
(with an infall parameter of f3(z = 3) = 0.38l 'Jg). Although 
the constraints on the mass density are relatively weak, we 
see that the constraints on the infall parameter remain con- 
sistent. Using these measurements to constrain the growth 
parameter, we find a value of f(z — 3) = 6/3 = 0.83±0.46. In 
addition we determine the combined parameter fas, which 
gives a measure of the growth parameter that is less depen- 
dent on the assumed underlying dark matter mass distri- 
bution. In this case we find a value of fas = 0.22 ± 0.12. 
These measurements are the highest redshift constraint on 
the growth parameter based on galaxy clustering analyses. 
We have shown that these measurements are consistent with 
the ACDM standard model, although given the uncertain- 
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ties on the measurements, they are also consistent with a 
number of other cosmologies. 

This work is one of the largest surveys of the galaxy 
mass distribution at z ~ 3 and paves the way for a number 
of lines of research, which will be followed in subsequent 
papers. In particular, the proximity of the data presented 
here to quasar sightlines will provide important constraints 
on the relationship between galaxies and the IGM at an 
epoch associated with significant interactions between the 
two. 
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